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Abstract. In 1960, Sobolev proved that for a finite reflection group G, a 
G-invariant cubature formula is of degree t if and only if it is exact for all 
G-invariant polynomials of degree at most t. In this paper, we find some ob- 
servations on invariant cubature formulas and Euclidean designs in connection 
with the Sobolev theorem. First, we give an alternative proof of theorems by 
Xu (1998) on necessary and sufficient conditions for the existence of cubature 
formulas with some strong symmetry. The new proof is shorter and simpler 
compared to the original one by Xu, and moreover gives a general interpreta- 
tion of the analytically- written conditions of Xu's theorems. Second, we extend 
a theorem by Neumaier and Seidel (1988) on Euclidean designs to invariant 
Euclidean designs, and thereby classify tight Euclidean designs obtained from 
unions of the orbits of the corner vectors. This result generalizes a theorem 
of Bajnok (2007) which classifies tight Euclidean designs invariant under the 
Weyl group of type B to other finite refiection groups. 



A main problem of numerical integration is to approximate the integral 



Here x is an n-dimensional coordinate vector and fi is a, probability measure on 
a domain Vl in M". We search for an approximation formula by taking a positive 
linear combination of the function values of / at specified points xi, • • • ,XM^ that 
is, 



We call (jl.ip a cubature formula. The values Wi are the weights and xi are the 
points of a cubature formula. To each formula we assign the set of functions for 
which it is exact. Most often this set is the space of all polynomials of degree no 
more than i; in this case a cubature formula is said to be of degree t. We refer the 
readers to the comprehensive monograph [T^l US] for the basic theory of cubature 
formula. 

A fundamental objective is to construct cubature formulas of large degrees with 
few points. The requirement that a given cubature formula is exact for polynomials 
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up to a certain degree can be reduced to the problem of solving a system of algebraic 
equations. In general, the larger the number of points or the degree of a cubature 
formula is, the greater the size of this system is. Sobolev [331 gave a celebrated 
criterion to diminish the size of the system to be solved. Namely, he proved that an 
invariant cubature formula is of degree t if and only if it is exact for all polynomials 
of degree at most t invariant under the group. This is known as the Sobolev 
theorem. The Sobolev theorem is widely accepted by the cubature community in 
analysis and related areas; for instance see [^ . 

Independent of the line of research in analysis and related areas, Goethals and 
Seidel 14, Theorem 3.12] developed the invariant theory of Chebyshev-type cuba- 
ture formulas on the sphere or spherical designs. As a generalization of spherical 
designs Neumaier and Seidel [SQ] considered cubature formulas on several concen- 
tric spheres called Euclidean designs. Bajnok [2 classified tight Euclidean designs 
whose points are the union of the orbits of the corner vectors of the group i?„, 
and in particular, he obtained several new tight designs. Here a Euclidean design 
is tight if it is minimal with respect to a lower bound for the number of points. 
To obtain the results, Bajnok [21 Proposition 14] essentially used the idea of the 
Sobolev theorem for i?„-invariant Euclidean designs, though he did not offer the 
name of Sobolev. It seems that some researchers in combinatorics and related areas 
do not fully recognized the Sobolev theorem [5]. 

In this paper we find some observations on invariant cubature formulas in con- 
nection with the Sobolev theorem. In Section 2 we explain Sobolev's invariant 
theory in detail. We also explain some basic facts related to Euclidean designs, 
e.g., a theorem of Neumaier and Seidel }20j which is well known in algebra and 
combinatorics. In Section 3 we give an alternative proof of famous theorems by 
Xu [28l Theorem 1.1, Theorem 1.2] on necessary and sufficient conditions for the 
existence of cubature formulas with radial symmetry. The original proof by Xu 
requires some tedious calculations and technical tools in numerical analysis like, 
Gaussian-Lobatto quadrature, Gaussian-Radau quadrature. Eventually it is long, 
and researchers in other areas may not be familiar with his proof. Whereas, our 
new proof is short and simple compared to the original proof. Moreover it gives a 
general interpretation of the analytically-written conditions of Xu's theorems, and 
so will be readable and acceptable for researchers not only in analysis, but also in 
other areas like algebra and combinatorics. In Section 4 we extend the theorem 
of Neumaier and Seidel to invariant Euclidean designs, and thereby classify tight 
Euclidean designs obtained from unions of the orbits of the corner vectors. This 
classification generalizes the result of Bajnok for other finite reflection groups. 

2. Preliminaries 

Let Hom;(R") be the linear space of all real homogeneous polynomials of total de- 
gree I in n variables. Let ViiW) = ELo Hom,(M"), P;(K") = X;l=o^ Homi_2j(M"). 
We denote by Harmi(K") the subspace of ^((R") of harmonic homogeneous poly- 
nomials of degree I. Let T'; ( A) , P;* ( A) be the space of functions which are the 
restrictions of the corresponding polynomials to A C R". 

Let G be a finite subgroup of the orthogonal group in R" and / £ 7';(R"). We 
consider the action of ct € G on / as follows: 
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A polynomial / is said to be G-invariant if it satisfies that 



We denote by Harmi(R")'^ the set of G-invariant polynomials in 'Pi(R"), 

Harmi(K."') respectively. 

A cubature formula (jl.ll) is said to be invariant under G, or G-invariant if the 
domain Q and measure /i of the integral are invariant under G and the set of points 
is the union of G-orbits and to each point of the same orbit an equal weight is 
assigned. The following is known as the Sobolev theorem. 

Theorem 2.1 ([25J). With the above set up, a G-invariant cubature formula is of 
degree t if and only if it is exact for every polynomial f G Pt(M")'^. 

The Sobolev theorem is widely accepted by the cubature community in analysis and 
related areas: In particular Russian mathematicians in analysis have developed the 
Sobolev theorem and employed it to construct many cubature formulas; for instance 
see [m m] . Xu [35] presented beautiful theorems on the existence and structure 
of cubature formulas for radially symmetric integrals. In section 3 we review his 
theorems in detail and give a short proof using the Sobolev theorem. 

Next let us explain a combinatorial object called Euclidean design which was in- 
troduced by Neumaier and Seidel [20] . Let AT be a finite set in R" . Let ri , r2 , • • • , 
be the norms of the vectors in X. For i = 1,2,- ■■ ,p we denote by S"^^ the 
sphere of radius centered at the origin, namely, S"^"^ — {x € R" | ||a;|| ~ r.i}, 
and let A^ = A n Sf ^^. The set A is said to be supported by p concentric 
spheres S = [J^^i S"^^. To each Si we assign the surface measure pi. Let 
IS*-'"^! = Jg„-i dp^{x), with the convention that Jg^^-^ f{x)dpi{x) = /(O) 

if5r' = {o}. 

Definition 2.2. With the same notations as in the above paragraph, we say A is 
a Euclidean t-design supported by S if there exists a positive weight function w{x) 
on X such that 



for every polynomial / e Vt{S). 

We can regard a Euclidean design as a cubature formula on some concentric spheres. 
Conversely a cubature formula for a class of integral with some symmetry is a 
Euclidean design (cf. [15l Lemma 3.1]). The following theorem by Neumaier and 
Seidel is well known in algebra and combinatorics. 

Theorem 2.3 (|20|). With the same notations as in Definition \2.'A the following 
are equivalent: 

(i) X is a Euclidean t-design with a weight function w. 



(ii) 'Exexwi^)fix) = for every f G ||a;|pJHarm/(M") with I < I < t,0 < 



In Section 4 we give a stronger theorem than Theorem 12.31 for invariant Euclidean 
designs, especially for researchers in algebra and combinatorics. 

Define p' = p — es, where £s = 1 if € S*, and es = otherwise. The dimensions 
of ViiS) and V'liS) are well known. 



af^f, 



V(T e G. 
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Theorem 2.4 ([III HI)- Let S C 

(1) d:iuyVi{s) = 



s5 + E-f'„-Mn!:n */^>v, 



dim Vi (M" ) = [''f) ifl<2p'-l. 

+ Eto' ("+!r-r') ^ *s I > 2p', 

(2) dim 7^(5) = <; Eto' ("+'r-r') »^ o'^'^: ^ > V, 

dimTTiM") = Elio ("+!r-r') < V - 1- 

The following lower bounds are known as the Fisher-type inequality for the size 
of a Euclidean design [6l [11] [ITl [18] ; the latter one is also called the MoUer bound. 

Theorem 2.5. (1) Let X he a Euclidean 2e-design supported by S. Then, 

\X\ > dimVeiS). 

(2) Let X be a Euclidean (2e — l)-design supported by S . Then, 



\X\ > 



2 dim7'*_i(S') - 1 if e is odd and e X, 

2 dim V*_ i{S) otherwise. 



A Euclidean ^-design is said to be tight if the equality holds in one of the bounds 
in Theorem 12.51 

Hereafter we assume G is a finite irreducible reflection group in R" . It is known 
that finite irreducible refiection groups are classified completely [9j. Let integers 
1 — mi < 7712 < • • ■ < f^n be the exponents of G (see PI Ch.V, §6 ]). 

Theorem 2.6 (|14)). Let G be a finite irreducible refiection group. Let qi be the 
dimension of }iaj:mi(W^)'~^ . Then, 



n ^ 



oo 

i=a i=2 



Note that for any x e M", the orbit x'^ is a spherical 7rt2-design in 14J. 

Let ai,a2, . ■ . tOLti be the fundamental roots of a reflection group G. The corner 
vectors vi,V2, ■ . ■ ,Vn are deflned by Vi _L a j if and only if i ^ j. We may assume 
= 1. We consider the set 



X{G,J)^ U 



where J C {1, 2, . . . , 77} and rk > 0. Let R denote the set of r^. 

Bajnok [2] found new tight Euclidean designs in X{Bn,J). In Section 4, using 
the theory of G-invariant harmonic polynomials, we extend the method of Bajnok 
to other reflection groups G, and classify the tight Euclidean designs obtained from 
X{G, J). 



3. XU'S THEOREM 

Let rj = {{x,y) G I < ^/x^ + y^ < 00}. Let W be a nonnegative weight 
function on [0, cxd) with flnite moments. We consider the radial weight function de- 
fined by W^(A/a?^ + y^) on fi. The following integral is said to be radially symmetric 
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or circularly symmetric: 



I[f]= / fix,y)W{y/x^+y^)dxdy 
Jn 

/•oo / /'27r 



i) (i) ^(^'^°'^^''''^^"^)^^)^^Wd^ 



To generalize a famous theorem by Verlinden and Cools [27 on the existence of 
cubature formula for radially symmetric integral, Xu proved the following theorem: 

Theorem 3.1. ([2H1 Theorem 2.1, Theorem 2.2]). (i) Let 

" ' (2j + a,)7r ^ .^^ { 2j + gQ^ 

2to + 2 



m + 1 -^-^ V 2m + 2 



(3.1) l2m+i(/) -Ao/(0,0) + 

TT ^ {2j + a,)TT . (2j+CT07r 

r- > Ai > / f i cos — : — , r j sm —— ; — 

i—l j—0 ^ 

where takes the value if m + i is even and the value 1 if m + i is odd. Then, 
In forms a cubature formula of degree 2n — 1 for I if and only if the following two 
types of conditions are satisfied: 

.00 L"/2J 

(3.2) / r'''+^W{r)dr = V X^r^ , j = 1, • • • , n - 1, 

L"/2J 

(3.3) ^ V2^(-1)^ = 0, J- L(^j + 3)/2j,... ,n-l. 

i=l 

(a) With the same symbol ai as in (i), let 

L(»+2)/2j 2m . /o - , ^ /o - , ^ \ 

T (f\ 27r ^ \ f (2j+cr, 7r . (2j+cr,)7r\ 

/ = — — > A^2^/ r,cos ,r^sm , 

m + 1 ■'^ ■'^ V 2m + 1 2to + 1 / 

i—l j—0 

where n = 2m — 1 or 2m. Then, I„ forms a cubature formula of degree 2n fori if 
and only if the following two types of conditions are satisfied: 

.00 L("+2)/2j 

(3.4) / r2j+iM^(r)dr = ^ \rl\ j = l, 
•^0 i=i 

L("+2)/2j 

(3.5) ^ A.rf(-l)'-0, J- L(n+1)/2J,--- ,n-l. 

i=l 



The aim of this section is to give a proof of Theorem 13.11 different from the 
original one by Xu: 

Proof of Theorem \3.1\ Let Di be the dihedral group of order 2i and f{xi,X2) be 
a Df-invariant polynomial. Using the polar coordinate system, we let f{xi,X2) = 
f{r cos 6, r sind). It is shown that / can be represented as a polynomial in two 
variables r^,r^ cos £9. We now consider the case where n = 2m in (i); the reader 
will easily see that the same argument as below works in the remaining cases. The 
conditions (|3.2p . p.3p respectively mean to substitute the monomials r^, • • • , r^"^^ 
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and r"+2 cos (n + 2)6*, r"+'* cos {n + 2)6, • • • , r^"-^ ^os (n + 2)9 into ([SII]). Thus the 
result follows by Theorem 12.11 □ 

Remark 3.2. (i) Xu's original proof of Theorem 13.11 is about 2 pages long only 
in the case of (i) with n — 2m. With this in mind, we tried a short proof using 
the Sobolev theorem. An advantage of our proof is the simplicity. Namely, the 
Sobolev theorem is the only advanced tool used in our proof, whereas, Xu's proof 
requires some technical and advanced tools in numerical analysis like, Gaussian- 
Lobatto quadrature, Gaussian-Radau quadrature. The proof by Xu also requires 
some tedious calculations. Eventually our proof is shorter and simpler than the 
original one. Another advantage of our proof: The conditions p.2|) . ()3.3p (or p.4p . 
(|3.5|) ) are considerably analytic, and so researchers in other fields like combinatorics 
and algebra will not be familiar with them. Whereas, our new proof gives a gen- 
eral interpretation of the above analytic conditions, and will promise researchers 
in these areas to enjoy Theorem 13.11 well. The authors hope that researchers in 
many different fields know the importance of Theorem 13.11 and will be more inter- 
ested in classical theories of cubature formulas developed in numerical analysis, (ii) 
Bajnok [l] Theorem 9] found a tight Euclidean t-design of which has the same 
structure of points as Xu's formula, as a generalization of a tight 4-design by Bannai 
and Bannai [4]. To do this, he implicitly used the same idea as in Theorem l3.11 for 
instance Eq. (10) in his paper corresponds to Eq. p.3|) (or Eq. p.5|) ) of our paper. 



4. Orbits of a reflection group as Euclidean designs 

In this section we classify the tight Euclidean designs obtained from X{G, J) for 
a finite irreducible refiection group G. A finite set X C M'' is said to be antipodal if 
X = —X. A tight Euclidean 2e-design has a weight function which is constant on 
each Xi [4], and so does an antipodal tight Euclidean (2e— l)-design [8]. Throughout 
this section we assume a weight function is constant on each G-orbit. 

First, we look at a stronger theorem than Theorem l2.3l for G-invariant Euclidean 
designs. A Euclidean t-design X is said to be G-invariant if X is a union of G-orbits 
and to each point of the same orbit an equal weight is assigned. 

Lemma 4.1. Let G be a subgroup ofO{W^). Let f be a G-invariant polynomial 
and x'^ be a G-orbit. Then, f{y) — f{z) for any y, z ^ . 

Proof. Straightforward. □ 
Let |G| be the order of a group G. 

Theorem 4.2. Let G be a subgroup o/0(R"). Let X — U^f^j^rfcX^, where Xk G 
5""^^ and rj, > 0. The following are equivalent: 

(1) X is a G-invariant Euclidean t-design. 

(2) Y.xex'^{^)M\^^ V{x) = for any e Harm,(R")'^ with I < I < t,0 < 

j<L¥J- 

Proof. For / e Harm;(]R''), the polynomial 

' ' sec 
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is an element of Harm/ (M'*)'^. Let w{x) — Wk for every x G x'^ . By Lemma I4TT1 for 
any / e Harm;(K'^), we have 

M 

J2 w{x)\\x\\'^^ f{x) ^J2^krl^ /(^z^^) 
xex k=i xt^xf 



fe=l ' ' geG 



4) 



x^\(p{xk) 



fc=i 

M 

^Wkrl^ ^ Vivkx) 

k=l xGx^ 

xex 



The result thus follows by Theorem 12.31 □ 

Remark 4.3. The radii rk are not necessarily mutually distinct in Theorem 14.21 
Goethals and Seidel [14] stated Theorem 12.11 for spherical designs. Theorem 14.21 
with all Tfc = 1 means the theorem of G-invariant spherical designs. The approach 
of using the orbits under subgroups of 0{M.'^) has long been considered [5]. The- 
orem 14.21 reduces the computational cost to check the strength of a G-invariant 
Euclidean design less than using Theorem 



From now on, let G be a finite irreducible reflection group, and Vk be a corner 
vector. Put = The following is immediate by Lemma 14.11 and Theorem 

K2\ 



Corollary 4.4. X{G, J) is a Euclidean t-design if and only if there exist Wk > 
and rk > such that the equation '^keJ '^k'f'^''^^ Nkf(vk) — holds for any f £ 
Harmi(R'')'^, where i,j are positive integers satisfying 1 < 2j + i < t. 



The dimension of Harmi(IR'*)"-^ is clear by Theorem 12.61 We can determine the 
basis of Harmi(M'')'-^ by harmonic polynomials / satisfying f{x^) = f{x) for each 
generator g of G. By the basis of Harmi(IR'')'-^ and CoroUarv 14. 4[ we know a 
necessary and sufficient condition for X{G, J) to be a Euclidean i-design. Bajnok 
[2] found an explicit such condition for the group Bn by using CoroUarv 14.41 For 
other groups, it is possible to give the conditions, but the statements are not simple. 
Therefore we do not write them in the present paper. 

Now, let us classify the tight Euclidean design obtained from X{G, J). For each 
group, we determine the possible maximum strength of X{G, J) for any J and 
radii R. Since the cardinality of is easily calculated, we can check whether the 
total size of a union of several orbits attains the fisher type inequality. For the set 
attaining the bound, we give its maximum strength by CoroUarv 14.41 

Hereafter let G M" be the row vector whose i-th entry is 1 and other entries 
are 0. Let Sn be the symmetric group. Define 

sym(/) := — ^ ^ 
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where (5„)/ := {g & Sn \ f{xO) = f{x)}. 



4.1. Group An- 
Dynkin diagram 



ai a2 as 
• • •— 



Exponents 

l,2,...,n 

Fundamental roots 

cti := Ci — em-i for 1 < i < n — 1. an ■= [a,a, . . . ,a,b] where a = (—1 + + l)/n 
and b = {n — 1 + \/n + 1) /n. 
Corner Vectors 

Vk = [cfe, . . . ,Ck, dk, ■ ■ ■ , dk] whose first k coordinates are equal to Cfe, and last n — k 
coordinates are equal to dk, where 

n+l-k + ^/n+l , -k 



Ck 



^k{n+l - fc)(n + 2 + 2Vn + 1) ^k{n+l - k){n + 2 + 2^Jn+l) 



Reflection group 

The reflection group A„ C 0(M") is generated by the following: 

r{ai) = [*ei, • • • *ei_i, *ei+i, *ei, *ei+2, • • • , *e„] for 1 < i < n - 1, 



r{an) 



I„-i-a?{^jj) -ab^j 

~abj 1 - 52 

where /„ is the identity matrix of size n, and j is the all-ones row vector. 
Orbits 

Let fi be the set of all vectors such that k coordinates arc equal to Ck, and other 
n — k coordinates are equal to dk- Let U2 be the set of all vectors such that fc — 1 
coordinates are equal to c^, and other n+1 — k coordinates are equal to d'^, where 

, n + 1 — k , —k — \/n + 1 



^k{n+l - k){n + 2 + 2y/n + T) ^k{n+l - k){n + 2 + 2^/nTl) 

Then the orbit v^" = t/i U t/2. Furthermore, we have Nk = ("J^) and u^f" = 

^n+l-k- 

Harmonic Molien series 

, ( l+t^ + ■■■ , if n = 2, 

= < l+t^ + + + ■■■ , if n = 3, 



G-invariant harmonic polynomials 

1. Degree 3. 

Note that dim Harms (R")^" = 1 for any n > 2. Harms (R")^" is spanned by the 
following: 

(i) n = 2. 

fi = x\ — 3xiX2 — 3a;ia;2 + x^. 

(ii) n = 3. 

3 3 
/s = sym(xf ) - -sym(a;ia;2) - -sym(a;ia;2a;3). 
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(ii) n > 4. 

J- / 3\ 3 , 2\ , 6(2 - \Jn + 1) , , 

/s = syni(a;f) - -^^—^sym{xxx^) + _ _ 3) ^y™(^i^2^3)- 

2. Degree 4. 

Note that dimHarm4(IR2)^2 = and dimHarm4(IR")^" = 1 for any n > 3. The 
following are ^n-invariant harmonic polynomials: 

/i4,i = sym(a;t) ^—'SYm.{x\x^^, 

n — I 

'14,2 = sym(.xiX2a;3X4), 

6 

/i4,3 = sym(xia;2) -sym(xia;2a;3). 

Harm4(R")"^" is spanned by the following: 
(i) n = 3. 

J4 = n.4,1 — Y3^'^" 



(ii) n > 4. 



where 



/4 = /i4,i + 04,2/14,2 + 04,3/14,3, 



_ 24(n + 2)(n^ - 5n - 12 + 4VnTT) 
""^'^ ~ (n - l)(n - 2)(n3 - 2n2 - 15n - 16) ' 

4(n + 2)(n2 - 2n - 7 - (n - l)\ArTT) 



C4,3 



(n- l)(n3 - 2n2 - 15n - 16) 



3. Degree 5. 

Note that dim Harms (M")'^" = forn = 2,3, and dim Harms (M")^" = 1 for any 
n > 4. The following are Sn-invariant harmonic polynomials: 



/i5,i = sym(a;^) - --^sym(x^a;^) + _ ^^^^ _ 2) sym(a;ia:^a;^), 

10 5 
/is,2 = sym(a;f ) -sym(a;ia;2) + -sym(a;ia;2), 

n — 1 n — i 

h5,3 = sym(a;ia;2a;3) ^sym(a;iX2a;3X4), 

n — 6 

/is,4 = sym(a;iX2X3a;4Xs). 

Harm4(]R")'^" is spanned by the following: 

(i) n = 4. 

, , , 17-20^5, , 10(18 + V5), 

J5 = "5,1 H "5,2 H 'l5,3- 

08 87 

(ii) n > 5. 

/s = /l5,l + 05,2/15,2 + 05,3/15,3 + 05,4/15,4, 
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where 

C5,2 
C5,3 = 



C5,4 



2n^ + 5n? - 21n - 90 - n(n + 6)y/nTT 
4n3 + 3n2 - 60n - 180 ' 
20(2n3 + 6n2 - 32n - 168 + (n^ - 8n + 12)\/nTT) 
(n - l)(n - 2)(4n3 + 3n^ - 60n - 180) ' 
120(n + 6)(n2 - lln - 78 + (2n2 - 2n + 12)VnTl) 



(n - l)(n - 2)(n - 3){An^ + Sn^ - 60n - 180) 

Substitute for G-invariant hetrmonic polynomials 

1. Degree 3. 
For n = 2,n > 4, 

A; — 

where 
For n = 3, 



2(n3 + 3n2 - 12n - 16 + (Sn^ - 4n - 16)Vrr+T) 
(n - l)(n - 3)(n + 2 + 2x/nTT) i 



2. Degree 4. 
For n > 3, 



where 



faivi) = ™, h{v2) = 0, /3(t;3) = 



= fc(n+l-fc) 



6(n + l)(n5 + 7n4 - 24n3 - 160n2 - 256n - 128 + 4(71^ - 20n2 - 48n - 32)s/n+l) 
{n - l){n - 2)(n3 - 2n'^ - 15n - 16)(n + 2 + 2^^7+1)2 ' 



6 ' ^=^ + 6 • 



3. Degree 5. 
For n > 4, 



/sl^^fe) = 77; ^ —3 05(n), 



(fc(n + 1 - A:)) 
where 

, , , 24(;( + 1)(2;)'' + + 50;/-' - 448;r' - 2144;)- - 3200;/ - 1536) 

05(n) = - 

(n - l)(n - 2)(n - 3)(4n3 + Zn^ - 60n - 180)(n + 2 + 2\/nTT)3 

24(n + l)(lln5 + SOn^ - 9671^ - 1120n2 - 2432n - 1536)\/^r+T 



+ 



(n - l)(n - 2)(n - 3)(4n3 + 3n2 - 60n - 180)(n + 2 + 2V^r+T)t 



, n + 1 V3(n+l)(2n-3) ^, n + 1 ^^(^ + l)(2n - 3) 
6 ' ^=^ + 6 • 

Theorem 4.5. There is no choice of J, R and w for which (A'(A„, J),w) is a 
Euclidean 6-design. 
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Proof. The polynomial of degree 6 

f{xi,X2, ...,Xn)= sym(a;ia;2) - ysym(x^^a;f ) 
is harmonic for any n> 2. We can calculate 

•^(^) = giin,k){k{k-n - l)(4fc2 - 4(n + + + 5n + 4).g2H + .gsH), 

where 

n{n +1) fn — \^ 



^) + 1 - fc)3(n + 2 + 2V^m)3 \k - 1 

g2{n) = 5(n^ + lln + 12 + (6n + 12)V?i + 1) 



53(n) = (n + \f{n + 2)(2n2 + 28n + 30 + (ISn + 30)V'^rfT). 
Note gi{n, fc) > for 1 < fc < n. Define 

F(fc) = k{k-n - l)(4fc^ - 4(n + l)k + n'^ + 5n + 4:)g2{n) + gsin). 
For a fixed n, we prove F{k) < for 1 < fc < n. We have 



dk 
where 



F{k) = 16 f fc - ^ ) (fc - - /3")52H, 



„ n + 1 v/2(« + l)(^-2) ^„ n + 1 v/2(n + l)(7i - 2) 

a = , p ~ h — . 

2 4 ' ^ 2 4 

If F(l) = F{n) < 0, F{{n + l)/2) < 0, and F{a") = F{l3") < 0, then F{k) < for 

all 1 < fc < n. Indeed for n > 2, 

F(l) = 

= -(n - l)(3n* + 27n^ + lOn^ + 26n + 60 + {15n^ + 15n^ + 60)\^TT) < 0, 
F{a") = F{(3") 

= -^n{n + lf(5n^ + 63n^ + 68n - 16 + (SOn^ + 60n)Vn+ 1) < 0. 
16 

F (^-y-) = - + l)^(7fi^ + 59n + 60 + (30n + 60)y/n+l) < 0. 

Therefore J^ke.J J^xgv^^ Wkrkf{x) < for any J, R and w. □ 

k 

Theorem 4.6. X(An,J) is not a tight Euclidean t-design except for the sets in 
Table 1. 

Proof. We prove only the classification of tight Euclidean 4-designs on two concen- 
tric spheres obtained from X(An, J). The other cases can be proved by a similar 
way. 

faivk) and 03 (n) are defined as above. Since 4>3{n) > for n > 1, fsivk) = if 
and only if n = 1 (mod 2) and k = (n-|- 1)/2. Clearly fsivk) > for A: < (n-|- l)/2, 
and fsivk) < for k > {n + l)/2. Therefore J must contain fci and fc2 such that 
fci < (n + l)/2 < /C2 by Corollary SH 

The size of a tight Euclidean 4-design on two concentric spheres is (n-|-l)(rt-|-2)/2. 
By noting that {n + l)(n + 2)/2 < A^3 = Nn-2 for ti > 5, we can determine 
J — {l,n — 1}, (or equivalently J — {2, n}) for any n > 2, or J = {1, 2} for n — 2. 
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\R\ 


t 


n 


J 






1 


2 


any 


{1} 


n = 1 


Wl = 1 




2 


any 


in} 


r,i = 1 


Wn — 1 




3 


3 


{2} 


r2 - 1 


102 = 1 




5 


2 


{1,2} 


ri = r2 = 1 


Wl ^ W2 = I 




5 


7 


{2,6} 


r2 = re = 1 


W2 = We = 1 


2 


4 


2 


{1,2} 


= 1, '^2 1 


Wl = 1, W2 = A 

'2 




4 


4 


{1,3} 


'^1 = 1' '^3 = 75 


Wl = 1, W3 = 27 




4 


4 


{2,4} 


''4 = 1, ^2 = 


7«4 = 1, 7«2 = 27 




4 


5 


{1,4} 


ri = 1, 7-4 = y/f 


Wl = 1, W4 ^ ^ 




4 


5 


{2,5} 


^■5 = 1, r2 = yi 


W5 = 1, W2 ^ I 




4 


6 


{1,5} 


ri = 1, rs = Vl5 


Wl = 1, W5 = -y 




4 


6 


{2,6} 


f6 = 1, ?'2 = Vl5 


W6 = 1, W2 = -gj 




5 


3 


{1,2,3} 


''I = ?'3 = 1, ''2 7^ 1 


Wl = W3 = 1, W2 = -Jt 




5 


5 


{1,3,5} 


ri = = 1, rs ^ 1 


Wl = W5 = 1, W3 = 



Table 1. Tight Euclidean t-designs from X{An, J) 



For n — 2, we can obtain tight Euchdean 4-designs on two concentric spheres as in 
Table 1. 

f^ivk), (t>4in), a, and /3 are defined as above. Note that ip4{n) 7^ 0, a > f , and 
j3 < n for any integer n > 2. Therefore l<a<n— l</3, (or equivalently 
a < 2 < P < n) holds by Corollary 14.41 The integers satisfying the condition are 
only n = 4, 5, 6. For n = 4, 5, 6, we can obtain tight Euclidean 4-designs on two 
concentric spheres as in Table 1. □ 



4.2. Group Bn. 
Dynkin diagram 

Oti a2 as O-n-l an 

• • • • • • ■ ■ 

4 

Exponents 
l,3,...,2n- 1 
Fundamental roots 

ai :— Ci — e^+i for 1 < i < n — 1 and a„ :— \p2eri- 
Corner Vectors 

= [1/Vfc, . . . , l/Vk, 0, . . . , 0], where Vk has k coordinates equal to l/Vk. 
Reflection group 

The reflection group Bn C 0(M") is generated by the following: 

r{ai) = [*ei, • • • *ej_i, *ei+i, 'e^, *ej+2, • • • , *e„] for 1 < i < n - 1, 

r(a„) = [*ei, • • • *e„_i, -*e„] . 

Orbits 

The orbit vj^" is the set of vectors with exactly k nonzero coordinates equal to 
±l/Vk. Note that uf" is antipodal and Nk = 2''{l). 
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\R\ 


t 


n 


J 




Wi 


1 


3 


any 


{1} 


n = 1 


Wl = 1 




3 


2 


{2} 


r2 = 1 


W2 = 1 




7 


2 


{1,2} 


ri = r2 = 1 


Wl = W2 — 1 


2 


5 


2 


{1,2} 


'"1 = 1, ?'2 7^ 1 


Wl ^ 1, W2 = pc 
Wi = l,W3 = ^^ 




5 


3 


{1,3} 


ri = 1, r3 ^ 1 




7 


4 


{1,2,4} 


ri = r4 = 1, r2 ^ 1 


Wl = Wi ^ I, W2 = ^ 


3 


7 


3 


{1,2,3} 


ri-1, ra-^g^, (r, ^rj) 


Wl = 1,W2^ ^,W3 = 



Table 2. Tight Euclidean t-designs from X{Bn, J) 



HEirmonic Molien series 

1 ^ f l + t'^ + t^ + --- , if n = 2, 

(1 - _ t6) . . . (1 _ i2n) \ l+t'^+t^ + ... , if n > 3. 

G-invariant harmonic polynomials 

1. Degree 4. 

Note that dini(Harm4(M")^") = 1 for any n > 2. The following is a i3„-invariant 
harmonic polynomial of degree 4: 

/4 = sym(a;i) ^sym(a;ia;2). 

n — 1 

2. Degree 6. 

Note that dim(Harm6(R2)^^) = and dim(Harm6(K")^") = 1 for any n > 3. The 
following is a i?„-invariant harmonic polynomials of degree 6: 

/e = sym(a;f ) - -J^sym{xlx^) + j--—^^—^sym{xlxlxl). 

Substitute Vk for G-invariant harmonic polynomials 

1. Degree 4. 



n — 1 



, , , 1 / /fc- 1 (k- 2) 
/6(t'fe) = 7^ 1-15 - + 30^ ^ 



2. Degree 6. 



fc2 y 71-1 (n- l)(n-2) y ■ 

Theorem 4.7 ( 2 ). There is no choice of R, J, and w for which {X{Bn, J),w) is 
a Euclidean 8- design. 

Theorem 4.8 {^). X(BmJ) is not a tight Euclidean t-design except for the sets 
in Table 2. 

Remark 4.9. We can also prove Theorems 14.71 and 14.81 bv the i?„-invariant har- 
monic polynomials. 
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4.3. Group Dn- 
Dynkin diagram 



ai 



a2 



as 



an-2 




Exponents 
1,3, ...,2n- 3,n- 1 
Fundamental roots 

a, := Cj — ej+i for 1 < z < n — 1 and a„ := e„_i + e„. 
Corner Vectors 

Vk = [1 /Vk, . . . , 1/ \/fc, 0, . . . , 0] , where Vk lias fc coordinates equal to 1/ for 1 < 
k < n-2. Vn-i = [l/v^, ■ • ■ , I/aAj, -1/v^] B.ndvn = [l/v^, • • ■ , 

Reflection group 

The reflection group Dn C 0(IR") is generated by the following: 

r{ai) = [*ei,---*ej_i,*ei+i,*ej,*ej+2,--- ,*e„] for 1 < i < n - 1, 

^(an) = [*ei, • • • *e„_2, -*e„, -*e„_i] . 

Orbits 

For 1 < fc < n — 2, w^" = u^". The orbit (rcsp. v^"i) consists of the 
vectors {il/y^}" with an even (resp. odd) number of negative coordinates. Note 
that v^" = —v^^i for odd n, and both u^" and v^-i are antipodal for even n. 
Furthermore, |A^„_i| = \Nn\ = 2""^. 
Harmonic Molien series 



(1 - _ f6) . . . (1 _ f2„-2)(l _ 



1 + 2*4 + t'^ + 3<8 + • • • , 

1 + t"* + + + 2*^ H 

1 + t"^ + 2t^ + 2t^ -\ , 

l + t^ + t'^ +f + 2t^--- , 
1 + + + 2t^--- , 



if n = 4, 
if n = 5, 
if n = 6, 
if n = 7, 
if n > 8. 



G-invariant harmonic polynomials 

1. Degree 4. 

Note that dim(Harm4(K4)^*) = 2 and dim(Harm4(M")^") = 1 for any n > 5. The 
following are Z)„-invariant harmonic polynomials of degree 4: 



/4 = sym(a;i) - ^^—^sym{xlxl). 

The following is a £>4-invariant harmonic polynomial of degree 4, which is linearly 
independent of f^: 

fi,2 = X1X2XSX4. 

2. Degree 5. 

Note that dim(Harm5(M^)^^) = 1 and dim (Harms (R")^-) = for any n ^ 5. The 
following is a Ds-invariant harmonic polynomial of degree 5: 

/5 = X1X2X3X4X5. 

3. Degree 6. 

Note that dim(Harm6(M*^)^<') = 2 and dim(Harm6(K")^'') = 1 for any n 7^ 6. The 
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following is a D„-invariant harmonic polynomial of degree 6: 

/e = sym(xf ) - ^^^sym(a;?a;^) + j^^—^^^^-—^sym{xlxlxl). 

The following is a Dg-invariant harmonic polynomial of degree 6, which is linearly 
independent of /e: 

/6,2(wfe) = a;ia;2a;3a;4a;5a;6. 

Substitute Vh for G-invariant heirmonic polynomials 

1. Degree 4. 

For 1 < fc < n - 2, 

/4K) = Kl-3'^-^ 

For k = n — l,n, 



k \ n—1 



fi{vk) = --■ 
n 



For n = 4, 



U,2{vi) = 0, fi,2{v2) = 0, /4,2(V3) = fiA'"^) = TZ' 

lb Id 

2. Degree 5. 
For n = 5, 

h{vi)=Q, h{v2) = Q, h{vz)=Q, h{vi) = --^, h{vi) = 



25\/5' 25\/5" 



3. Degree 6. 

For 1< A: < n - 2, 



/6(^^fc) = 72 1-15 7 + 30^ ' 



k^ \ n-1 (n- l)(n-2) 

For k = n — l,n, 

For n = 6, 

/6,2(-(^l) = 0, /6,2('i^2) = 0, /6,2(W3)=0, /6,2(-(^4) = 0, 

/6,2('f5) = -2T6' feMve) = 2T6- 

Theorem 4.10. There is no choice of J, R and w for which {X{Dn, J),w) is a 
Euclidean S- design. 

Proof. The following is a £>„-invariant harmonic polynomial of degree 8: 

/g = sym(a;i) ^sym(a;ia;2) + — ^sym(a;ia;2). 

n—1 n—1 



For 1 < A: < n - 2, 
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li?l 


t 


n 


J 






1 


7 


8 


{2,7} 


r2 -- 


= rr = l 


W2 




= 1 




7 


8 


{2,8} 


r2 -- 


= rg = 1 


W2 


= W8 


= 1 


2 


5 


6 


{1,5} 


n = 


1, r5 1 


Wi = 


1, W5 


_ 9 
8rf 




5 


6 


{1,6} 


ri = 


1, re 7^ 1 


Wi = 


1, ^6 


_ 9 
8rJ 




7 


4 


{1,2,3,4} 


ri =r3 = 


'"4 = 1, ?'2 7^ 1 


Wi =W3 = 


i(;4 = 


1, W2 = 4 

'^5 



Table 3. Tight Euclidean t-designs from X{Dn,J), where 
n or n — 1 G J 



and for k — n — l,n, 

Therefore fsivk) > for all fc. □ 

Theorem 4.11. Assume J contains n or n — 1. Then X{Dn,J) is not a tight 
Euclidean design except for the sets in Table 3. 

Proof. By Theorem 14 . 1 01 and £)„-invariant harmonic polynomials, a proof is similar 
to that of Theorem SH □ 

Remark 4.12. The tight Euclidean designs in Tables 1,2,3 are already known in 

[illlSllTllHllTn] 

Remark 4.13. For each G = F4, H3, H4, Eq, Ej, Eg, by checking the cardinality 
of a union of several v^, we can prove X{G,J) is not a tight Euclidean design 
except for known tight spherical designs [2TJ [23] . 

5. Concluding remarks 

In this paper we found some observations on invariant cubature formulas and 
Euclidean designs in connection with the Sobolev theorem. First, we gave an alter- 
native proof of celebrated theorems by Xu on necessary and sufficient conditions 
for the existence of cubature formulas with radial symmetry. The new proof is 
much shorter and simpler compared to the original one by Xu. Thus researchers 
in analysis will realize again the importance of the Sobolev theorem. Moreover 
our proof gives a general interpretation of the analytically-written conditions of 
Xu's theorems, and so will promise researchers in algebra and combinatorics to be 
more familiar with Xu's theorems. Second we extended the Neumaier-Seidel theo- 
rem to invariant Euclidean designs, and thereby classified tight Euclidean designs 
obtained from unions of the orbits of the corner vectors. The classification gener- 
alizes Bajnok's theorem to other finite reflection groups beside groups of type B. 
Bajnok's theorem and results obtained in Section 4 may imply that invariant cuba- 
ture formulas of high degree could hardly exist. Xu [22] pointed out, however, that 
the general Lie groups has been used for studying cubature formulas in a different 
setting - cubature rules on the fundamental domain of the group, which are for 
exponential or trigonometric functions - and they yield Gaussian type cubature for 
algebraic polynomials of very high orders; for instance see [16], [22] for details. We 
believe this direction of research in analysis will also motivate the study of cubature 
formulas in other areas of mathematics. 
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